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Abstract 

We set up a general scheme to resum ladder diagrams for the quark- 
anti-quark potential in A/* = 4 super- Yang-Mills theory, and do explicit 
calculations at the next-to-leading order. The results perfectly agree with 
string theory in AdS^ x when continued to strong coupling, in spite of 
a potential order-of-limits problem. 
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1 Introduction 



The potential between static charges in gauge theories is defined through 
the expectation value of a rectangular Wilson loop. In the AdS/CFT con- 
text, Wilson loops are directly related to string worldsheets and at strong 
coupling their expectation values obey the area law [H El The static po- 
tential in A/" = 4 super- Yang-Mills theory (SYM) can thus be computed at 
strong coupling by evaluating the area of a minimal surface in AdS^ x 
with appropriate boundary conditions \T\ [2] . This calculation has led to one 
of the first quantitative predictions of the AdS/CFT duality flj. Recently 
much progress has been achieved towards a non-perturbative description of 
the static potential at any coupling. The main tools are diagram resumma- 
tions [4j and quantum integrability of the AdS/CFT system [5j. The former 
approach is less general and is applicable only in a corner of the parameter 
space, but is technically simpler and can be derived from first principles. 
Apart from giving another insight into the inner working of the AdS/CFT 
duality, this approach can be helpful in finding direct gauge-theory interpre- 
tation of various quantities appearing in TBA (Y- functions, driving terms in 
the TBA equations etc.). 

The approach of |4] is based on the resummation of ladder diagrams. The 
Bethe-Salpeter equation that resums ladders for rectangular Wilson loop was 
derived in [6j and predicts the \/A scaling of the static potential, consistent 
with the string calculation at strong coupling. The limit that emphasizes 
ladder diagrams involves analytic continuation in the scalar coupling of the 
Wilson loop [4], and once the same limit is taken on the two sides of the 
duality, the results completely agree [4] , in spite of a potential order-of-limits 
problem. Identification of the ladder approximation with the leading order 
of an expansion in a small parameter [4J makes it possible to systematically 
compute correction^ Our aim will be to set up a general scheme of doing 
ladder resummation order by order. We will explicitly compute the leading 
correction to the simple ladder approximation. In particular we will check if 
the order-of-limits problem start to aflFect the comparison to string theory at 
strong coupling, once the NLO corrections are taken into account. 

^ See [7 for an early discussion of how to improve the ladder approximation. 
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2 Coulomb potential 



The Wilson loop operator in A/" = 4 SYM theory is defined as 



^^(C) =trPexp 



^ ds {Afj,x^ + i^in^\x\) 
c 



(2.1) 



It describes an infinitely heavy W-boson associated with the symmetry break- 
ing by a Higgs vev in the direction of the six-dimensional unit vector n: 
($/) = vni. Rectangular T x L loop with the scalar couplings iii and n2 on 
the two sides of the rectangle defines the Coulomb potential between W and 
W with the symmetry-breaking vevs 111^2: 

ln{W{CT.L,n,,n,))-^^^J^. (2.2) 

The Coulomb charge a{9^X) is a function of the 't Hooft coupling A = Qyu^ 
and the angle between the symmetry-breaking vevs: cos 9 = iii • n2. 
To the leading order in perturbation theory, 

a{9, A) = ^ (1 + cos 9) + (A^) , (2.3) 

where the first term in the brackets comes from the gluon exchange and the 
second term from exchanging scalars. The limit, which emphasizes the ladder 
diagrams, requires an analytic continuation in 9: 

9 = i^. (2.4) 

When is real and large, the scalar exchange becomes dominant. In the 
scaling limit [4], 

i9^oo, Ae^- fixed, (2.5) 

the scalar contribution is of order one, while the gluon contribution remains 
small because we should simultaneously take A ^ 0. It is easy to see that any 
scalar line connecting the two sides of the Wilson loop is enhanced by a factor 
of ni-n2 = COS 9 ^ e^, while adding a gluon line or an internal vertex produces 
an uncompensated coupling A. The scalar ladders (diagrams without internal 
vertice^ are thus the only diagrams that survive in the limit i9 ^ oo [4j. 

^ There are also diagrams without internal vertices, where propagators connect points 
on one and the same line. However such diagrams cancel, since a separate straight Wilson 
line is supersymmetric. 
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Resummation of the ladder diagrams can be reduced to solving a simple 
Schrodinger equation [6j. 

It is convenient to introduce a rescaled coupling 

A = - (cos + 1), (2.6) 
2 



which is of order one in the limit (2.5), (2.4), and to consider the Coulomb 
charge as a function of A and A. The ladder expansion can be regarded as 
series in e"^, or equivalently in A: 

oo 

a(A,A) = ^a„(A)A- (2.7) 

Our goal is to compute ai. To make the presentation self-contained, we will 
also review the leading order of the ladder approximation that computes ao. 



3 Ladder diagrams 

3.1 Bethe-Salpeter equation: general structure 

To calculate the n-th order correction to the Coulomb charge, in principle 
it is enough to identify all diagrams that have order O(A^) in the scaling 



limit (2.5). But the number of such diagrams is infinite, and a better way to 
organize the calculation is to first compute 2PI diagrams and then generate all 
other diagrams with the help of the Bethe-Salpeter equation. The number of 
2PI diagrams at each order is finite. In addition, the Bethe-Salpeter equation 
automatically exponentiates the result, making it easy to take the logarithm 



in (2.2). 



The kernel of the Bethe-Salpeter equation K{s' ^s^t' ^t) is the sum of 2PI 
diagrams whose end-points lie within the intervals {s^s') and {t^V) on the 
two anti-parallel sides of the Wilson loop. The sum of all diagrams, not 
necessarily 2PI, with end-points within the intervals (0, S) and (0,T), which 
we denote by r(5', T), then satisfies the Bethe-Salpeter equation (fig. [T]): 

S s' T t' 

r(5,r) = l + J ds' J ds J df J dtK{s\s]t\t)r{s,t). (3.1) 
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= 1 + t r(s,t) C^s\s^rxi 



Figure 1: The Bethe-Salpeter equation. 



The Coulomb charge is related to the large-distance behavior of r(S', T): in 
the units where 1/ = 1, 

(3.2) 

The equation can be simplified by differentiating with respect to S and T: 

S T 



f ds J dtK{S,s]T,t)r{s,t). (3.3) 



dSdT 



It is convenient to introduce new variables (fig. [T]): 

y^S + T, x^S-T, u^S-s, v^T-t. (3.4) 

Owing to the translation invariance, the kernel does not depend on the center- 
of-mass coordinate y: 

K{S,s]T,t) = K{u,v]x), (3.5) 
and the Bethe-Salpeter equation becomes 

y+x y-x 
Q2Y Q2Y r r 

— ^ - — ^ = du dv K(u^v]x)V(y -u-v^x -u + v). (3.6) 
ai/ ox^ J J 

^ 

The equation should be supplemented with the boundary condition 



T{y,x) 



= 1. (3.7) 

y^\x\ 



According to (3.2), the Bethe-Salpeter wavefunction T{y^x) grows expo- 
nentially with at least when x = 0. Under the assumption that at small but 
finite X the dependence on y is still exponential, we arrive at the following 
ansatz: 

T(y,x)^ijix)e^ , (3.8) 

where a is the Coulomb charge that we want to calculate. This ansatz can be 
valid only at |x| « as is does not carry any information on the boundary 



condition (3.7). At the leading ladder approximation it is possible to make a 



more rigorous argument and derive the behavior (3.8) from the first principles 
[8j. We will not go into such a detail here, and will simply assume that (3.8) 
holds for sufficiently large y. This assumption will be sufficient for finding 
both a and ip{x). 



Substituting the ansatz (3.8) into the Bethe-Salpeter equation (3.6), and 
sending y to infinity while keeping x finite, we get an integro-differential 
equation for 



oo oo 



dx'^ 





J du J dv e '^^''^''^K{u,v]x)ijj{x-u + v) ^^ip. (3.9) 



The Coulomb charge a is the smallest eigenvalue of this linear problem. 



3.2 Summing up ladders 

The Bethe-Salpeter kernel can be computed in perturbation theory: 

oo 

K{u, v]x)^Y. ^""^niu, V] x), (3.10) 

where each order contains a finite number of 2PI diagrams. Each is a 
function of A. 

The only 2PI diagram of order O(A^) is the scalar exchange: 

\5{u)5{v) 

Ko{u,v]x) = — — . (3.11) 
47r2(x2 + l) 



With this approximation, the eigenvalue problem ( |3.9 ) reduces to an ordinary 
Schrodinger equation: 
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which has been analyzed in great detail in [6l [HI [9l Hj. In particular, the 
eigenvalue of the Schrodinger operator can be systematically analyzed at 
large and at small A. For completeness, we review the results below. 

The weak-coupling expansion of ao is not exactly a power series: the 
Taylor coefficients contain logarithms. Such a non-analytic behavior arises 
because of the IR divergences of perturbation theory for the static potential 
[10]. The logarithms first show up at two loops and can be resummed 
using a version of renormalization group equation [TT]. Both logarithmic 
and finite terms can be calculated from the Schrodinger equation. To the 
three-loop accuracy j^: 



A A2 Ae7-i A3 / Ae7 Ae^-i 7 7r^\ ^^^^^ 
q/q(A) = — + — - In + In In + • • • , (3.13) 



where 7 is the Euler constant. 

At large A, the eigenvalue of the Schrodinger problem (3.12) is given by 
the classical energy at the minimum of the potential 

ao = — -! + .... (3.14) 

TT 

The second term is the energy of the zero-point fiuct nations. The strong- 
coupling asymptotics of the ladder diagrams perfectly agrees with the clas- 
sical calculation in string theory, as soon as the correct range of parameters 
is identified f4j. The angle 9 defines the boundary conditions for the string 
worldsheet on in the AdS^ x geometry. Once the angle is analytically 
continued to complex values and i9 is taken to infinity, the area of the mini- 
mal surface in AdS^ x grows exponentially, as e^/^, reproducing the ladder 



result (3.14). It would be very interesting to compare the second term with 
the one-loop correction due to string fiuctuations, which can be calculated 
along the lines of [12]^ 

The wavefunction in the strong-coupling limit is highly peaked near zero: 

1 

^o(x)-(^)'e-4f^ (A -00). (3.15) 

The ladder diagrams are thus getting densely packed at large A [7J. In- 
deed, the average (x^) ^ A~^/^ determines the mean square distance between 



^We thank Nadav Drukker for the discussion of this point. 
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(a) 



(b) 



Figure 2: Comparison of the exact eigenvalue of the Schrodinger equation (3.12) 



(sohd blue line) with the three-loop weak-coupling result ( |3.13 ) (dashed magenta 
line) and one-loop strong-coupling approximation ( 3.14[ ) (dot-dashed yellow line). 



consecutive rungs of the ladder. At strong coupling this distance becomes 
parametrically small. 

In fig. [2] the numerical solution for in the full range of A is compared 
to the weak and strong coupling approximations above. 



4 Beyond ladders 

To compute the NLO contribution to the Coulomb charge, we first need to 
calculate the 0(A) correction to the kernel - the sum of 2PI diagrams pro- 
portional to A, A^cos^ and A^cos^^ (there are no 2PI diagrams of higher 
order in A). The relevant diagrams are shown in fig. [sj The resulting cor- 



rections to the Schrodinger equation (3.12) are non-local, which definitely 



complicates the problem. On the other hand, we do not need to solve the 
equation exactly, but can apply ordinary quantum-mechanical perturbation 
theory. As we shall see, the problem really simplifies when A is very large 
(at strong coupling) and when A is very small (at weak coupling). 

Although the rectangular Wilson loop is not a supersymmetric observ- 
able, the diagrams in fig. [3j when combined together, feature spectacular 
cancellations. In the Feynman gauge the gluon and the scalar propagators 
are equal. The diagram (a) then renormalizes the effective coupling, and it 
has been already taken into account by including the 1 in the definition of 
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vwwwwv 



(a) (b) (c) (d) (e) 

Figure 3: Bethe-Salpeter kernel at the next-to-leading order. 



p-k 



k-q 

ivwvwvl 



p-q 
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(a) (b) 
Figure 4: The H-diagram. 



A, eq. (2.6). The sum of the remaining three diagrams can be considerably 



simplified. 

We start with the H-graph, fig.|3](e). The numerator in the corresponding 
momentum integral (the momentum assignment is shown in fig. [4^), can be 
written as 

-{k + q){2p-k-q)^ -2/ + k'^ + q^ + {p- kf + {p- qf - {k - qf. (4.1) 

This is the standard Passarino-Veltman reduction. All the terms but one 
contain a vanishing propagator. Those terms combine into 



vwvww 



+ 




The last term exactly cancels with the X-diagram in fig. |3]l. 

The Y-graph in fig. [3](c) can be reduced to the scalar diagram with the 
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derivatives along the contour, acting on the external lines: 



(4.3) 



after which the external leg of the diagram can be integrated by parts. The 
fact that the Y-diagrams gives a total derivative has been observed many 
times in similar calculations [H [13]. Upon integration by parts only the 
boundary terms survive: 




(4.4) 

The last two diagrams cancel off the four middle terms in (4.2). As shown in 



the appendix |A| the bubble diagram cancels with the self-energy corrections 
in fig [3|3. We are thus left with a single contribution, given by the second 
derivative of the scalar H-graph (fig. lib): 



(4.5) 



Ki{u^ V] x) = Ai7(xi, . . . , X4), 
where the Laplacian is defined as 



A = (ai + a2)'. 



(4.6) 



and 



(4.7) 

The diagram is taken in the following kinematics: 

xi = (u, 0, 0, 1), X2 = (0, 0, 0, 1), X3 = (x, 0, 0, 0), X4 = (x + V, 0, 0, 0). (4.8) 

We do not know if one can express the H-diagrams in terms of known 
functions. This is a difficult double-triangle two-loop integral, for which we 
were unable to find any closed analytic expression. The best we could do 
is a double-integral representation, presented in appendix [Bl However, for 
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analyzing the strong and weak coupling limits it will suffice to know the 
asymptotic expressions of the H-diagram, which are relatively easy to get 
from the Feynman integral directly. 

The problem of finding the Coulomb charge at the NLO boils down to 
solving the linear problem (3.9) with the kernel given by the sum of (3.11) 
and (4.5). Doing this exactly is rather difficult, as the Schrodinger operator 
becomes integro-differential, and the eigenvalue a enters the equation non- 
linear ly. However, A is small and we can apply the usual formulas of quantum- 
mechanical perturbation theory. But as we do not know the leading-order 
wavefunctions explicitly, we will restrict ourselves to the two limiting cases 
of small A and large A. 



5 Strong coupling 



At large A the problem significantly simplifies, because then a is large, and 



the integral transform in (3.9) is highly peaked at small u and because 



of the exponential suppression factor. In particular, the scale of variation of 
the wavefunction is proportional to A"^/"^, cf. (3.15), and is thus much larger 



than the scale at which the integrand gets exponentially suppressed, which 
is ~ A~^/^. We can thus neglect the shift in the argument oi hj v - u 
and reduce the problem to the Schrodinger equation with a local effective 
potential 



(5.1) 



Moreover, the kernel itself changes on the distances of order one, and we thus 
need to know its expansion at small u and v. This corresponds to the regime 
when the end-points of the H-diagram collide pairwise (fig. [T]). The relevant 
limit is computed in eq. (B.8), with logarithmic accuracy: 



2567r6(x2 + l)^ 



(5.2) 



Substituting this result into (5.1), we get for the effective potential: 



A 



AA^ In a 



+ 



47r2 + 1) 27r^a^ (x^ + 1)' 



(5.3) 
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At strong coupling, when the eigenvalue is given just by the depth of the 
potential well, we get: 



a= 1 ^ . 5.4 

^ V 2A ' 



The second term in the brackets is the correction to the ladder approximation. 



Since we have not determined normalization of the logs in (5.2), the answer 
is calculated with the logarithmic accuracy. It should be possible to more 
accurately determine all the constants and calculate the non-logarithmic term 
of order A/A, but we will not attempt to do it here. 

It is a non-trivial fact that the corrected Coulomb charge, as a function 
of A and is still proportional to A^/^ and thus can be directly compared to 
the string calculation at strong coupling. Provided that In A can be traded 
for in the strong-coupling limit, we get: 

o^.^ef (l-2i9e-^). (5.5) 

5.1 Comparison to string theory 

In string theory, the Wilson loop vev is given by the area of the minimal 
surface bounded by anti-parallel lines on the boundary of AdS^^ and inter- 
polating between two points on with angular separation 9. The solution 
for the minimal surface was found in [Ij. The area can be written in the 
parametric forn][^ 

„ = ^■y-^l^-(^)'<f (5,0) 

tt/cv 1 - fc^ 

^ = 2\/2F - 1 (5.7) 

where K = K{k^) and E - E(k'^) are the complete elliptic integrals of the 
first and second kind. The angle 9 is real when the elliptic modulus satisfies 
k'^ < 1/2. Analytic continuation to imaginary 9 corresponds to real fc^ > 1/2. 
We are interested in the limit i9 ^ oo, which corresponds to taking fc^ ^ 1. 



^The same result can be extracted from the minimal surface bounded by a cusp j3l[T2]. 
We follow the notations and conventions of [12J. 
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Expanding (|5.6|), (|5.7|) in 1 - fc^ with the logarithmic accuracy, we get: 



a - -^(l.l^ln(l-fc^).O(l^ (5.8) 

^ = -ln(l-P) + lnl6+^-^i^^ ln(l-P) + 0(l-P). (5.9) 
Consequently, 

a^^e^l-2^e-^). (5.10) 

in an agreement with ( |5.5[ ). 

As emphasized in [4J , taking A large and then sending i9 to infinity is not 
quite the same as keeping A small and sending A to infinity. We see however 
that potential order-of-limits ambiguities do not affect the agreement of the 
ladder resummation with string theory even at the NLO order. 



6 Weak coupling 



The first order of perturbation theory in A for the eigenvalue equation (3.9) 



produces the following answer for the correction to the Coulomb charge: 

2A ^ ^ ^ 
5a ^ — dx J du J dv e~^^'^^'"'^ilJo{x) ipo^x - u + v) Ki(u^v^x)^ (6.1) 



where Ki is first correction to the kernel from (4.5)-(4.8). At small A the 
unperturbed wavefunction is given by 



(A-0). 



(6.2) 



Since ao, given by eq. (3.13^, is parametrically small in A, all exponential 



to one, which leads to 



factors in (|6.l|) decay very slowly and to the first approximation can be set 
!ads to 

oo oo oo 

5a = \}? J dx J du J rf'^ Ai7(xi, . . . , X4), (6.3) 



where Xi are expressed through v and x according to (4.8), or according to 
fig. [1] We got the expected result, that to the leading order in A it is enough 
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Figure 5: The dangerous region of integration in the H-diagram. 



to integrate the H-diagram along the Wilson loop. There are no ladders 
to sum, as extra rungs of the ladder are accompanied by extra powers of A. 



However, as we shall demonstrate in a moment, the integrals in (6.3) actually 



diverge. To understand how the divergence is cut off, we need to go one step 



back to the all-loop expression (3.9). The exponential suppression factors in 



the exact formula becomes important on the scale 



t 







A- 



(6.4) 



This scale imposes an effective IR cutoff on the integrals in (6.3). 



The integral (6.3) is convergent by power counting, which indicates that 
nothing special happens when all points go to infinity simultaneously. The 
divergence comes from the configuration in which X2 and X3 are a finite 
distance from one another, while Xi, X4 go to infinity (fig. [5]), and from a 



symmetric configuration in which Xi, X4 stay close, while X2, X3 go to infin- 



ity. In this regime we can use the long-distance approximation (B.9) for the 



H-diagram. Moreover, we can apply the Laplacian only to the logarithmic 



pref actor, since the derivatives acting on the other terms in (B.9) raise pow- 
ers in the denominator, and produce IR convergent integrals. Then, using 
Alnx23 ^ 4/^23 = 4/(x^ + 1), we can easily integrate over x: 



-t-cx 

/ 



dx 
+ 1 



= TT. 



(6.5) 
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We are thus left with the v and z integration: 

^ ^ oo oo 

r r r 

5a = -2 — -— du dv — = Tprj ^ • (6.6) 

4%M J J z^[z^ + (u-z^f][z^ + {v-z^f] 

The factor of two takes into account the identical contribution from xu 1- 
Doing the u and v integration first, we arrive at 

r d'^z (tt 

da = -- — — / — — — + arctan-- . 6.7 
(27r)7 J z^T? \2 |z|/ ^ ^ 



The integral is most easily computed by changing variables to p = v and 
(/9 = arctan(z^/|z|): 



The logarithmic divergence at small p is fake. It arises because we have 
initially assumed that p » 1, and neglected, for example, the separation 
between the two sides of the Wilson loop. The neglected contributions reg- 
ularize the divergence at some p ~ 1. On the contrary the IR divergence 
at large p is real, and is regularized only by the resummation of an infinite 
number of ladder diagrams. The resulting cutoflF, up to a numerical factor, is 



given by (6.4). Taking this relationship into account, we get the final answer 



for the correction to the Coulomb charge (3.13): 



<^« = ^AA2logA (6.9) 



Taken at face value, the result (6.9) is a three-loop single-logarithmic 



correction, which combines with the single-log terms in (3.13). There is also a 
finite (9(AA^) term, that we have not calculated. Clearly, the higher orders in 
the ladder expansion will produce single-log contributions of the form A^ A In A 
and A^ In A, which originate, for instance, from the purely gluonic H-diagram. 
It would be interesting to compute these log-enhanced contributions, which 
complete the perturbative three-loop calculation of the Coulomb charge in 
the NLLO. 
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7 Conclusions 



The leading-order ladder resummation can be done also for the Wilson loops 
on X [4J. In this configuration maps to a cusp with the opening 
angle (/p, where (f is the angular separation of the Wilson lines on the sphere 
[Hj. Our calculation corresponds to the limiting case ip ^ n. It would be 
interesting to evaluate the NLO correction for general 99, especially since the 
Schrodinger problem is exactly solvable at 99 = 0. This latter case corresponds 
to the straight Wilson line, with the scalar coupling jumping from to cos^ 
at some point on the contour. 

It would be also interesting to reproduce the ladder limit (2.5) from the 
TBA equations, that in principle describe the Coulomb charge at any value 
of 9 and A [5j. The TBA equations are highly non-linear, and in general are 
difficult to solve, but when 1} ^ 00 they should reduce to a simple Schrodinger 
equation (3.12). Such a dramatic simplification is not unconceivable, as 1} 
plays the role of a chemical potential in TBA. Large real chemical poten- 
tial (real 9 actually corresponds to an imaginary chemical potential) should 
strongly emphasize the degrees of freedom with positive charge, which per- 
haps can be taken into account by rearranging the ground state of the auxil- 
iary spin system. Whatever the underlying reason for the drastic reduction in 
complexity in comparison to the general case, in this corner of the parameter 
space the TBA equations should be solvable. 
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A Cancellation of the self-energy corrections 



Here we will show that the one-loop correction to the propagator (fig. 



exactly cancels with the bubble graph in the right-hand-side of eq. (4.4). Both 



of these diagrams diverges, and it is necessary to regularize them exactly in 
the same way, to make sure that no finite leftover survives the cancellation. 
We use the usual dimensional reduction scheme. 

In the dimensional reduction regularization, the one-loop scalar propaga- 
tor is given by [8]: 



r(f-i) 



AP(f-l) 



The bubble diagram then gives: 



(A.l) 



2Ar3(f-i) 



AP(f-l) 



^20-2(r-^)^-2 167r^ (L> - 3) (4 - r20- 



(A.2) 

which exactly cancels the second term in (A.l), even before taking the D ^ 4 
limit. 



B H-diagram 

The following auxiliary integral is useful: 



f (z^ + 



(z^ + /x^)^ + jy'^)^ {z - y + a) 



where uj± are the roots of the quadratic equation 



— , = f^(ln^^-ln^M, (B.l) 



(B.2) 



With the help of this formula, we can get for the H-diagram a double-integral 
representation over two sets of Feynman parameters: 



1 1 

H ^ / rf/^i / 

(47r)' J J 



(B.3) 
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where now uj^ are the roots of 

Cjf - [/^l (1 - + (1 - /^l) ^24 + ^1^2^?2 + (1 - ^1) (1 - ^2) ^34] ^± 

+ /^l/^2 (1 - /^l) (1 - /^2) X?4^23 = 0, (B.4) 

and 

a = /^iXi4 - /^2^23 - ^34. (B.5) 

It is also true that 

0? = /^l (1 - 1^2) ^?3 + /^2 (1 - /^l) ^24 + ^1^2^?2 + (1 - /^l) (1 - /^2) ^34 

-/^l (1 - /^i) - /^2 (1 - /^2) ^23- (B.6) 

B.l Short-distance asymptotics 

First we consider the case when |xi2| ~ IX34I « |x23|. The H-diagram is non- 
singular in the limit Xi2,X34 0, and on dimensional grounds we expect the 
leading term of the form const lx\^. As 1/^23 is annihilated by the Laplacian 
at ^2^X3, we will need the next term in the expansion in X12, X34. This term 
appears at the quadratic order and contains an extra logarithmic enhance- 
ment. The large logarithm ln(x^2/^23) ^^ses from the region of integration 
where both integration points simultaneously approach Xi^2 01* ^3,4- When 
both z and y are close to Xi, X2, each of the two "long" propagators can then 
be approximated by 1/^23, and we are left with the integral 

^ f ^ "^y"^ 2 2+(l'2-3,4). (B.7) 

The integral is the Fourier transform of equal to In 2:^2/2^71"^ ■ Taking 
into account the contribution from the other coincidence limit, we obtain: 

const ^n'^^^ + ^M^^^ 

^23 ^ /I ^23 

B.2 Long-distance asymptotics 

Next we consider the case when |x23| « \xi2\ ~ \xs4\ (fig. [5|). When X23 0, 



the ^-integration in (4.7) is log-divergent at ^ ^ X2,3. The divergence is cut 
off at ~ X23, and we can thus replace 

(z-xif {Z-X2Y {z-X'iY 



^^_liix|3 r 

4%8 J (z-x.? (z-xo? (z-x.?' ^ ' ' 
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